Experimental Investigation of an Eight-Qubit Unit Cell 
in a Superconducting Optimization Processor 
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A superconducting chip containing a regular array of flux qubits, tunable interqubit inductive cou- 
plers, an XY-addressable readout system, on-chip programmable magnetic memory, and a sparse 
network of analog control lines has been studied. The architecture of the chip and the infrastructure 
used to control it were designed to facilitate the implementation of an adiabatic quantum optimiza- 
tion algorithm. The performance of an eight-qubit unit cell on this chip has been characterized by 
measuring its success in solving a large set of random Ising spin glass problem instances as a function 
of temperature. The experimental data are consistent with the predictions of a quantum mechanical 
model of an eight-qubit system coupled to a thermal environment. These results highlight many of 
the key practical challenges that we have overcome and those that lie ahead in the quest to realize 
a functional large scale adiabatic quantum information processor. 



I. INTRODUCTION 

Of the known paradigms of quantum computing, those 
related to quantum annealing (QA)^"^ arc unique in that 
they leverage what ought to be a robust natural phe- 
nomenon - the tendency for physical systems to seek 
and remain in a low energy configuration. These meth- 
ods have received considerable theoretical attention of 
late with claims and counter-claims concerning their ul- 
timate utility^"-'^^. In contrast, they have received rela- 
tively little experimental attention. This is a significant 
concern as data are needed to help clarify the impact of 
environmental noise upon such approaches to quantum 
computing^^""'^^. Experiments have been performed us- 
ing nuclear magnetic resonance in molecules^°, but the 
prospects for scaling this approach to larger systems are 
limited. At least one proposal for an adiabatic quantum 
optimization (AQO) processor based upon superconduct- 
ing flux qubits exists^-'^'^^ and rudimentary experiments 
have been performed upon related small-scale devices 
with fixed qubit and coupler parameters^^^^^. Experi- 
mental investigations of QA in large-scale systems have 
only ever been performed upon solid state samples in 
which the experimenter has no control over the individual 
couplings between Ising spins'^^'^^. Consequently, there 
are many open questions regarding how to implement a 
large-scale programmable QA information processor in 
practice. 

The purpose of this article is threefold: First, we ad- 
dress some of the practical questions regarding how to 
design a scalable superconducting AQO information pro- 
cessor. The answers to these questions will then serve 
as a motivation for the architecture of the device that 
we have fabricated. Second, we present experimental 
results from a unit cell on one such chip consisting of 
eight flux qubits^®, sixteen in situ tunable inductive in- 
terqubit couplers^^, an XY-addressable high fidelity read- 
out architecture^*^, an array of in situ programmable fiux 



storage devices addressed by a single fiux quantum de- 
multiplexing circuit'^^, and a sparse network of analog 
control lines. The data demonstrate that the unit cell 
can be used as a computer for solving Ising spin glass 
problems. Third, we compare experimental data to the 
results of numerical simulations in order to highlight the 
fact that, when run very slowly with respect to the adia- 
batic limit, the performance of the unit cell is influenced 
by its tendency to thermalize to an environment. 

This article is organized as follows: Section II contains 
a mathematical description of how the physics of an Ising 
spin glass can be mapped onto superconducting hardware 
composed of rf-SQUID fiux qubits and couplers. Sec- 
tion III contains a brief review of the hardware that has 
been been fabricated and key calibration data from one 
of the eight-qubit unit cells. This particular unit cell was 
used to solve a large number of Ising spin glass prob- 
lem instances and the results have been summarized in 
Section IV. A dynamical model of the unit cell coupled 
to an environment in thermal equilibrium is introduced 
in Section V. Section VI contains a brief discussion of 
the key results of this article and presents a series of im- 
portant open questions that have been motivated by this 
work. The conclusions are summarized in Section VII. 
For convenience, we have provided a glossary of abbre- 
viations used in this article in Appendix A. A table of 
specific example problem instance settings referred to in 
Sections IV and V has been placed in Appendix B. 



II. MAPPING ADIABATIC QUANTUM 
OPTIMIZATION ONTO HARDWARE 

We begin by mapping a particular class of optimiza- 
tion problems onto a scalable processor architecture that 
uses superconducting fiux qubits. Significant emphasis 
will be placed upon minimizing the number of unique 
time- dependent control signals, thus making efficient use 
of the limited number of externally supplied biases that 
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can possibly be routed to such a processor. On the other 
hand, it will be assumed that unique time-independent 
control signals can be generated locally on chip us- 
ing a scalable form of programmable magnetic memory 
(PMM) that has been described in detail elsewhere^ 
Let there be an optimization problem of the form 
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where Si = ±1, —l<hi< +1, and —1 < Kij < +1. For 
any given set of hi and Kij there exists at least one opti- 
mal solution Sgs that minimizes the objective E. Finding 
Sgs for such a system of coupled variables can be NP- 
hard^^. Equation (1) can be mapped onto an Ising spin 
glass (ISG) Hamiltonian^'^ 
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where a (\ are Pauli matrices acting upon spin i and 
Eq is a convenient energy. When cast in this form, |sgs) 
represents the ground state of the ISG. The objective is 
to design a physical system that reliably finds |sgs). 

One possible means of finding |sgs) is via the adiabatic 
method described in Ref. 2. In an idealized picture of 
AQO, the Hamiltonian of the system at arbitrary time t 
can be expressed as that of a quantum Ising spin glass 
(QSG). In such a system, there are pairwise magnetic 
interactions between spins i and j oc Kij and each spin 
i is subjected to a local longitudinal magnetic field oc hi 
and a global transverse magnetic field oc T{t): 
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One finds the lowest energy solution to an optimization 
problem encoded in the quantities hi and Kij by guiding 
the physical system through an evolution path described 
by T{t) subject to the following constraints: 

r(0) > h„K,j , 

r{tf) < h„K,j , 

where tf > is the run time of the algorithm. Here, 
it is assumed that the physical system readily reaches 
its ground state jllo) at i = and, in the best of cir- 
cumstances, will remain in the ground state during the 
course of the evolution, thus yielding \sgs) &t t = t j . Un- 
derstanding the conditions under which adiabaticity is 
violated, thus yielding final states other than |sgs), is a 
matter of considerable debate at the moment^"^*^. 

It was recognized in Ref. 21 that one could implement 
a QSG using a network of rf SQUIDs in which some are 
designed to fulfill the role of qubits and others to fulfill 
the role of tunable inductive interqubit couplers. The 
physics of the two lowest lying states of an isolated flux 
qubit i can be captured by an effective Hamiltonian 
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tunneling energy, respectively. Here, |/^| represents 
the magnitude of the supercurrent flowing about the rf- 
SQUID loop, $f is the external flux bias applied to the 
qubit loop about which the supercurrent flows and 
represents the qubit degeneracy point. Upon compar- 
ing terms with like symmetry in Hamiltonians (3) and 
(4), one can see that Aq/2 = i?o(i)r(t). In order to 
implement an AQO algorithm in which T{t) is altered, 
one must provide a means for tuning A^ in situ. This 
can be accomplished by incorporating at least one dc- 
SQUID loop into an rf-SQUID body^S-^^-s? since we 
have used the compound-compound Josephson junction 
(CCJJ) rf SQUID in our experiments^®, let the external 
flux bias applied to this device be represented by ^ccjj- 
This leads to the substitution A^ - 
tonian (4). Rather unfortunately, 

of in ^ny compound junction rf-SQUID flux qubit: 
' q I 
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, in Hamil- 
is also a function 



K|'(*^ccjj)|- Therefore, one does not naturally ob- 
tain orthogonal control of and Ag via $f and $ccjj' 
respectively. This drawback is by no means unique to 
the CCJJ design, as it occurs in all flux qubits with in 
situ tunable A^ that have been reported on in the litera- 
ture to date. This intrinsic lack of orthogonal control of 
the two defining parameters of the flux qubit, |/^| and 
Aq , has driven some of the design choices made in imple- 
menting the QSG discussed herein. 

Assuming that one has a set of flux qubits for which 
|^9(*^ccjj)| ^nd Ag($j?j,jj) are nominally identical, a gen- 
eral time-dependent Hamiltonian for a system of induc- 
tively coupled flux qubits can be written as 
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where the time dependence of the qubit parameters |/^| 
and Ag is implicitly driven by ^^^^^{t). In addition, there 
could be time-dependence in the individual qubit flux 
biases and in the interqubit effective mutual in- 

ductances Mij via their control biases iji^)'^^- 

One of the virtues of the AQO algorithm is its sim- 
plicity, which provides some significant advantages when 
considering how to implement a processor. Note that 
Hamiltonian (5) is a function of a global bias ^^^^^{t) to 
which all qubits are uniformly subjected. This is de- 
sirable for two reasons: First, Hamiltonian (3) specifies 
that the transverse field Aq{t)/2 ~ Eo{t)r{t) be uniform, 
which then corresponds to all Ag($^j.jj) being identical. 
As a corollary, this renders all | identical. Sec- 

ond, one can provide uniform •I'ccjj to multiple qubits 
using a single global current bias line, as opposed to one 
bias line per qubit. This scenario is depicted in Fig. 1, 
in which the current bias /ccjj {t) drives the fluxes ^ccjj i^ 
a multitude of qubits. In doing so, one substantially re- 
duces the number of external biases that must be applied 
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to the processor, thus significantly improving the scaling 
of this particular architecture. Such design choices are 
most likely necessary in order to realize multiqubit pro- 
cessors that contain more than lO's of qubits. 

Given the ideal Hamiltonian (3) and the flux qubit- 
based device Hamiltonian (5), one must now map the 
problem specification, embodied by hi and Kij in Eq. (1), 
onto external control parameters. This can be accom- 
plished by scaling Hamiltonian (5) and then performing 
a term-by-term comparison to Hamiltonian (3). A con- 
venient scaling factor is the interqubit coupling energy 
when a coupler is set to provide maximum antiferromag- 
netic (AFM) coupling, Mij {<^co,ij) = Mafm- Define this 

I 1 2 

energy scale as JafmC^) = -^afm ■ Rearranging 

Hamiltonian (5) yields 
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FIG. 1: (Color online) Mapping of ISG problems, specified 
by a set of values denoted as hi and Kij , onto superconduct- 
ing hardware. Two qubits, two \ compensators, and one 
interqubit coupler are shown. A global current bias /ccjj(f) 
provides the fluxes ^^^^^{1) that drive the annealing process 
to multiple CCJJ rf SQUID flux qubits. Interqubit coupling 
is mediated by tunable mutual inductances Mij oc Ktj that 
are controlled by static fluxes $co,ij. Qubit bodies subjected 
to the sum of static flux biases and time-dependent flux 
biases driven by a global current bias Ig{t). The latter signals 
are mediated to each qubit via tunable mutual inductances 
Ali oc hi that are controlled by static fluxes i. 
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In order to solve a particular optimization problem, 
hi and Kij must be time-independent. According to 
Eq. (6c), one must hold all Mij constant during oper- 
ation. This is convenient as it obviates the application of 
individually tailored time-dependent fiux bias signals to 
each interqubit coupler. Rather, one need only apply a 
static control signal to each coupler, as depicted in 

Fig. 1, that can be provided by PMM. On the other hand, 
according to Eq. (6b), one must apply time- dependent 
qubit flux biases of the form 



= *° + /^. X MAFM\int)\ 



to render hi time-independent. Thus, it is necessary to 
provide a custom-tailored time-dependent control signal 
plus a static offset to every qubit. The static compo- 
nent $° can be provided by PMM. As for the time- 
dependent component, providing these signals with one 
external bias per qubit would not constitute a scalable 
approach for building a multiqubit processor. Rather, 
one can take advantage of the fact that, according to 
Eq. (7), all qubits must receive a control signal with 
the same time-dependent shape but with custom-tailored 
time- independent scale factors (x hi. A scalable architec- 
ture for providing these signals is depicted in Fig. 1 and 
has been further expounded upon in Ref. 31. Here, a 



single global current bias Ig{t) — a\lP{t)\, where a is 
a convenient scale factor, is coupled to multiple qubits 
via in situ tunable mutual inductances of magnitude 
Mi = hiAlAFu/ct. The very same type of device that is 
used to provide in situ tunable interqubit coupling^^ can 
be retooled to provide coupling between flux qubits and 
a global bias line. Moreover, each Mi can be controlled 
with a static flux bias 3>f^ ^ provided by PMM. We will 
hereafter be refer to this architecture as persistent current 
{\lq\) compensation, as it is a means of compensating for 
changing |^q(i)| such that the ISG problem specified by 



(7) hi and Kij remains on target throughout annealing. 



To summarize up to this point, a prescription for im- 
plementing AQO to solve ISG problems using a net- 
work of inductively coupled CCJJ rf-SQUID flux qubits 
has been presented. A problem specified by a set of hi 
and Kij can be embedded in the hardware using time- 
independent interqubit couplings controlled by PMM and 
time-dependent qubit fiux and CCJJ biases. The qubit 
flux bias signals can be supplied using a combination of 
static flux offsets provided by PMM and a single global 
signal Ig {t) that is applied to each qubit through in situ 
tunable couplers that are also controlled by PMM. The 
CCJJ bias can, in principle, also be provided to all qubits 
simultaneously using a single global control signal. 



III. DEVICE ARCHITECTURE AND 
CALIBRATION 

With the mapping of the AQO algorithm onto the scal- 
able architecture completed, we now turn to a high level 
description of a superconducting chip whose architecture 
embodies that algorithm. All of the principal compo- 
nents of the processor, namely the qubits^*, couplers^^, 
readout'^", and PMM^^ have been described in detail in 
other publications. As such, we will only provide brief 
summaries of the important points as pertaining to the 
functioning of the collective system herein. 

As stated previously, we have incorporated CCJJ rf- 
SQUID flux qubits in our design^®. A schematic of an 
isolated flux qubit with the two external bias controls 
relevant for this study is shown in Fig. 2(a). This partic- 
ular qubit is robust against fabrication variations in the 
Josephson-j unction critical currents and facilitates the 
homogenization of the net critical current among a popu- 
lation of such qubits. This device also contains an induc- 
tance (L) tuner that can be used to compensate for vari- 
ations of qubit inductance due to fabrication and from 
tuning the interqubit couplers^^. To each qubit we have 
added an compensator, as introduced in Sect. II. We 
have provided PMM to flux bias the two minor lobes of 
each CCJJ^^, the L tuner, the compensator, and the 
qubit body ($^) for all qubits on the chip'^"'^. 

The chip that was used for this study was composed 
of sixteen eight-qubit unit cells that were tiled on a 4 x 4 
square grid. To limit the scope of this article, we focus 
upon a single unit cell near the center of the chip. A 
discussion of the complete processor, with multiple unit 
cells acting in concert, will be reserved for a future pub- 
lication. A schematic layout of the unit cell is shown in 
Fig. 2(b). The qubits are depicted as extended horizontal 
and vertical loops with 16 compound Josephson junction 
(CJJ) couplers located at the intersections of the qubits. 
Additional couplers at both cxtrema of the qubit bodies 
connect them to qubits in neighboring unit cells to the 
left, right, top and bottom. These latter couplers were 
set to zero coupling to isolate the single unit cell for this 
study. Each qubit was connected to its own quantum flux 
parametron (QFP) enabled readout^^'^°. Not shown in 
Fig. 2(b) are the PMM elements, the demultiplexing tree 
for addressing the PMM, or the analog bias lines used to 
calibrate and operate the unit cell. 

A simplified representation of the hardware connectiv- 
ity is depicted in Fig. 2(c). Here, qubits and couplers 
correspond to the vertices and edges of a graph'^®. As 
with the unit cell schematic, this eight-vertex graph can 
be tiled to the left, right, top, and bottom in order to 
generate larger graphs. Each vertex is connected to a 
minimum (maximum) of 4 (6) other vertices, depending 
upon the number and arrangement of unit cells used to 
form a larger graph. This hardware does not provide full 
connectivity within a population of N qubits in which 
each physical qubit is connected to TV — 1 other physi- 
cal qubits. This limitation can potentially be overcome. 




FIG. 2; (Color online) Device schematic, (a) A single CCJJ 
rf SQUID with the two time-dependent biases relevant for 
this study. The flux bias ^%^y^{t) drives the annealing pro- 
cess. The global -compensation bias is provided by Ig{t). 
(b) Schematic layout of the eight-qubit unit ceU. Qubits are 
shown in grey. Example readout (RO), CCJJ, L tuner (LT), 
j/q I compensator (IPC), internal coupler (ICO), and portions 
of external couplers (XCO) have been noted, (c) Graph rep- 
resentation of the hardware connectivity. Qubits are repre- 
sented by vertices (solid dots) and couplings by edges (solid 
lines). Black (red) edges correspond to ICO (XCO) couplings. 



at the cost of reducing the number of unique vertices 
< N, by using a fcrromagnctically coupled chain of phys- 
ical qubits to form a single logical qubit, as suggested in 
Rcf. 22. 

The chip used in these experiments was fabricated on 
an oxidized Si wafer with Nb/Al/Al2 03/Nb trilayer junc- 
tions and four Nb wiring layers separated by planarized 
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plasma enhanced chemical vapor deposited Si02 . A scan- 
ning electron micrograph of the fabrication cross-section 
can be found in Fig. 4(b) of Ref. 28. The Nb metal lay- 
ers are referred to as BASE, WIRA, WIRE and WIRC, 
from bottom to top, respectively. Flux qubit wiring was 
primarily located in WIRB and consisted of 2-/im-wide 
leads arranged as an approximately 900- /im- long differ- 
ential microstrip located 200 nm above a ground plane in 
WIRA. Coupler wiring was primarily located in WIRC, 
stacked on top of the qubit wiring to provide inductive 
coupling. PMM flux storage loops were implemented as 
stacked spirals of 13-20 turns of 0.25-/im-wide wiring with 
0.25 ^m separation in BASE and WIRA (WIRB). Stored 
flux was picked up by one-turn washers in WIRB (WIRA) 
and fed into transformers for flux-biasing devices. Ex- 
ternal control lines were mostly located in BASE and 
WIRA. Resistors that were used in the PMM demulti- 
plexing circuit were made from a TiPt layer referred to 
as RESI. All of these control elements resided below a 
ground plane in WIRC. The ground planes under the 
qubits and over the PMM/external control lines were 
electrically connected using extended vias in WIRB so 
as to form a nearly continuous superconducting shield 
between the analog devices on top and the bias circuitry 
below. Transformers for biasing qubits, couplers, QFPs 
and dc SQUIDs were enclosed in superconducting boxes 
with BASE and WIRC forming the top and bottom, re- 
spectively, and vertical walls formed by extended vias in 
WIRA and WIRB. Minimally sized openings were placed 
in the vertical walls through which the bias and target 
device wiring passed at opposing ends of each box. This 
design reduced most on-chip parasitic crosstalk to a neg- 
ligible level. 

An optical image of a unit cell completed through the 
processing of WIRB is shown in Fig. 3(a). One can dis- 
cern the trenches in which the qubit wiring resides, where 
the bottom is formed by ground plane in WIRA that 
is electrically connected to thick vertical walls that are 
formed by extended vias in WIRB. A number of PMM 
flux storage and pickup loops are visible as spirals and 
washers, respectively. The PMM circuitry is not visible 
in a completed chip as it resides below patches of ground 
plane in WIRC that are electrically connected at their 
perimeter to the top of the vertical walls defining the 
qubit trenches. Thus the qubits (PMM) reside above (be- 
low) the nearly contiguous contiguous shielding layer in 
a completed chip. Note that much of the coupler wiring 
is absent from this image as it resides in WIRC atop the 
qubit wiring in WIRB. 

Tiling of the eight-qubit unit cell to make a larger pro- 
cessor is explicitly demonstrated in Fig. 3(b). One can 
discern a 4 x 4 array of eight-qubit unit cells. The par- 
ticular unit cell used in this study was in the third col- 
umn and third row. Outside of the field of this image are 
four blocks of on-chip hysteretic dc SQUIDs, with pickup 
loops oriented in-plane and orthogonal to the wafer, that 
are used as magnetometers to measure local fields during 
active field compensation, as well as wiring channels that 




FIG. 3: (Color online) (a) Optical image of an eight-qubit unit 
cell completed through the processing of WIRB. Qubits, la- 
beled as go — >■ 97, reside within the trenches formed by ground 
plane in WIRA and extended vias in WIRB. PMM elements 
are visible as spirals and washers between qubits. (b) Opti- 
cal image of a 4 x 4 array of eight-qubit unit cells completed 
through the processing of WIRB. Unit cell shown in (a) is 
enclosed within the dashed box. 



run to wirebonding pads. There are 128 qubits, 352 cou- 
plers and 128 readouts on this processor. These analog 
components would require a total of 1632 unique biases 
to operate if no effort was made to develop a scalable 
control architecture. It would be impractical to build 
such a processor with so many independent external bi- 
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ases due to circuit density and wirebonding constraints. 
However, our complete processor, consisting of the analog 
components enumerated above and 992 PMM elements, 
requires only 84 pairs of differentially-driven external bi- 
ases to calibrate and operate. 

The chip was enclosed in a superconducting Al shield 
(Tc ~ 1.2 K) that was inside two coaxial /i-metal shields. 
All three shielding layers and the sample were thermal- 
ized to the mixing chamber of a pulse tube dilution re- 
frigerator. The refrigerator itself was located inside two 
room temperature coaxial /i-metal shields. The magnetic 
field in the vicinity of the chip was minimized at 3.9 K 
using triple-axis Cu compensation coils located outside 
of the Al shield. The array of on-chip dc-SQUID magne- 
tometers surrounding the qubit circuitry was used to infer 
the magnetic field vector near the processor. The current 
in the compensation coils was controlled via feedback 
implemented in software. Once the magnitude of the 
three-dimensional magnetic field vector had been mini- 
mized, the chip was thermally cycled to > 9.5 K to expel 
trapped flux from the Nb chip {Tc « 9.1 K) and then 
cooled to base temperature T w 19 mK. Thereafter, the 
current applied to the compensation coils was removed. 
The residual magnetic field was estimated to be '--^ 2 nT 
(3 nT) perpendicular (parallel) to the plane of the chip. 

External biases applied to the chip were provided by a 
custom-built 14-bit 128-channel differentially driven pro- 
grammable 50 MHz arbitrary wave form generator. The 
signals generated at room temperature were sent into 
the refrigerator using individually shielded twisted pair 
wires that were connected to a combination of lumped- 
element and copper powder filters at the mixing cham- 
ber. To minimize the effect of environmental noise and 
since the AQO algorithm does not require resonant exci- 
tation of any element of the circuit, we chose to restrict 
the bandwidth of these bias lines to 5 MHz. Signals were 
then routed onto Sn traces on a printed circuit board 
(Tc w 3.7 K) and passed onto the chip via Al wirebonds. 

A summary of the flux qubit parameters that are rele- 
vant for this study is presented in Fig. 4. Data have been 
plotted versus normalized flux, where $o = ^/2e is the 
flux quantum. The reader is referred to Ref. 28 for a de- 
tailed description of the CCJJ rf-SQUID flux qubit and 
the experimental methods used to obtain (*5'ccjj)|- Wc 
have obtained Aq{^^^--) using two methods. First, we 
used the single qubit Landau-Zcner (IQLZ) method^^ to 
yield 2 x lO'"^ Hz < Aq/h < 10^ Hz. The lower and upper 
bound were set by experimental run time constraints and 
the bandwidth of our external control lines, respectively. 
Second, we ran two-qubit Landau-Zener (2QLZ) experi- 
ments on pairs of qubits that were strongly coupled with 
Mij = — 2.82pH. The Hamiltonian for a pair of qubits, 
subject to zero net flux bias, can be expressed as 



1 r 

2 



+M,,\If\\ma^a<^K (8) 



Data from a 2QLZ experiment yields a spectral gap 
which can be compared to that obtained by solving for 
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FIG. 4: (Color online) Measured and predicted qubit param- 
eters as a function of ^ccjj- (a) Magnitude of the persistent 
current \lq\- (b) Tunneling energy Aq. Solid curves are pre- 
dictions from the CCJJ rf SQUID Hamiltonian using inde- 
pendently calibrated device parameters. 



the eigenenergies of Hamiltonian (8). The range of g 
that could be inferred from measurements was also lim- 
ited by the run time and bandwidth constraints cited 
above. In the first variant of 2QLZ, we locked the tar- 
get (f>^^jj-biases of the qubits together. Assuming that 
Ai = Aj = Aq, one can use Hamiltonian (6a) to solve 
for Aq(<i>^j.jj) given the independently calibrated 5($^j,jj), 
|/f and \lj (<I>^(,jj)|. This procedure proved ef- 

fective for extracting 5 x lO^'Hz < Aq/h < 2 x 10® Hz. 
In the second variant of 2QLZ, we targeted ^ccjj ~ 
— 0.6235 $0 for qubit i (a bias for which we have cali- 
brated A,; using the IQLZ method) and scanned the tar- 
get 'I'jfj.jj of qubit j. We then used Hamiltonian (6a) to 
solve for Aj ~ Aq($^^jj) given the independently cali- 
brated g(*^eii): A,(-0.6235 3>o), |/f (-0.6235 $o)|, and 



CCJJ / 

$|5'j.jj)|. This procedure proved effective for extract- 



ing 4 X 10® Hz < Aq/h < 3 X 10^ Hz. 

Prior to collecting the data shown in Fig. 4, we had 
fully exercised the CCJJ minor lobe and L-tuner biases 
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on each of these qubits in order to homogenize |/^(<I>^^jj) 



at «'ccjj/*o = -0.6146 and -0.6325 among the eight 



qubits. Residual discrepancies between |/^| and Aq of 
different qubits at arbitrary ^^^.jj were attributed to small 
variations in the mutual inductance between the global 
CCJJ current bias line and the individual qubits. 




0.4 0.5 

So,,j/*0 




FIG. 5: (Color online) (a) Mutual inductance Mij versus cou- 
pler flux bias ^ CO, ij for ah sixteen interqubit couplers in the 
unit cell. Solid curves are the theoretical Mij obtained using 
the mean of the individual coupler ht parameters, (b) Mutual 
inductance Mi versus [-compensator flux bias $7^,;. Solid 
curve is the theoretical Mi obtained using the mean of the 
individual |/g [-compensator flt parameters. 



The theoretical predictions from a quantum mechani- 
cal Hamiltonian describing the CCJJ rf SQUID, as given 
by Eqs. (5a)— 7-(5d) in Rcf. 28, arc shown as solid curves in 
Fig. 4. The model results were obtained using the means 
of the independently calibrated CCJJ rf SQUID param- 
eters, namely the body inductance ibody = 339 ± 3 pH, 
CCJJ loop inductance iccjj = 26 ± 1 pH, critical current 

= 3.355 ± 0.040 AiA and capacitance Cq = 182 ± 4fF. 



Given that the experimental data agree with a multilevel 
quantum mechanical of the CCJJ rf SQUID whose input 
parameters were independently calibrated, we are confi- 
dent in our identification of these devices as flux qubits^®. 



The interqubit couplers and [/|'[ compensators on this 
chip were calibrated according to the methods described 
in Ref. 29. The experimental data are summarized in 
Fig. 5. From the set of calibrations for all 352 interqubit 
couplers on the chip, we determined that the smallest 
maximum AFM coupling was Mafm = 1-88 ± 0.03 pH, 
which was readily achievable for all sixteen couplers 
within the unit cell used in this study. We chose this 
value of AfAFM to define the energy Jafm introduced in 
Sect. II. Likewise, from the complete set of 128 [/^[ com- 
pensators on this chip we determined that the maximum 
usable AFM coupling was max \Mi\ = 0.100 ± 0.002 pH. 
We found it convenient to use the scale factor a ~ 
2 X Mafm /max [Mi [ = 37.6 to scale the compensa- 
tion signals, as discussed in Sect. II. 
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FIG. 6: (Color online) Veriflcation of readout fidelity for the 
eight-qubit unit cell, (a) Example results from sequentially 
preparing the eight qubits into each of the 2* — 256 possible 
final spin configurations and recording the number of correct 
reads for 1000 repetitions, (b) Probability of observing a given 
number of errors within 1000 repetitions, as taken from the 
data in (a). Data have been fit to a Poisson distribution with 
mean number of errors A = 0.36. 



Finally, to characterize the readout fidelity of the unit 
cell, we set all interqubit couplers to zero coupling and 
then sequentially prepared and read all 2* = 256 pos- 
sible final spin configurations \s) — \sq S1...S7), with 
Si = ±1. This was accomplished by initializing the 
qubits in the presence of a bias vector of the form 
e (X ( So si ... S7 ). We prepared each configuration 
and read the state of the unit cell 1000 times using two 
repeated reads of each dc SQUID'^". The results of one it- 
eration of this fidelity check arc shown in Fig. 6(a). Here, 
we have adopted a shorthand notation by converting each 
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binary word \s) into a decimal number | State), where 
State ^ ^^2^-^ . 

i=0 

In the majority of instances the readout yielded the ex- 
pected result 1000 times. We never observed more than 
3 errors per 1000 repetitions. Multiple iterations of this 
check revealed that results for those states that appear to 
be less than perfect in Fig. 6(a), such as State=143, were 
not reproducible. Further analysis indicated that these 
infrequent errors were not due to problems in prepar- 
ing \s). Rather, the errors were most likely generated by 
the statistical nature of the dc SQUID switching times^". 
The probability of observing a given number of errors 
within 1000 repetitions, as obtained by taking a his- 
togram of the data in Fig. 6(a), is shown in Fig. 6(b). 
These results are well described by Poisson statistics with 
a mean number of errors A = 0.36 in 1000 repetitions. 
We defined the threshold for a number of counts to be 
statistically significant as three times the square root of 
the variance, 3^/A w 2. One can then deem the fidelity of 
our readout of the eiglit-qubit unit cell in these particular 
experiments to have been 1 — 3-\/A/1000 = 0.998. 



IV. EXPERIMENTAL PROCESSOR 
PERFORMANCE 

To test the performance of the eight-qubit unit cell, we 
generated a set of 800 unique random ISG problem in- 
stances whose hi and Kij were drawn from a distribution 
of 15 evenly spaced values within the set 

h,, K,j e [-1 - 6/7 . . . -t- 6/7 + 1] . (9) 

We then took advantage of the symmetry of the physical 
layout of the qubits and couplers shown in Fig. 2b and 
generated 8 permutations of each of the 800 progenitor 
problems by rotating and reflecting their hi and Kij em- 
beddings about the center of the unit cell. The result 
was a test set of 6400 problem instances, each requiring 
a unique embedding, that were then posed to the hard- 
ware. Note that while other researchers have chosen to 
focus upon particular classes of Ising spin problems, such 
as random 3-SAT^°, MAX-chque'^^ , and exact cover^'^^, 
we have found our random ISG problem set to be of great 
use as it encompasses a broad range of problem classes 
that could be posed to an AQO processor in practice. 

The PMM programming required to set up the proces- 
sor was rather involved and the details of its implementa- 
tion are beyond the scope of this articlc'^^. Nonetheless, 
it is worth recognizing that, for each problem instance, 
the interqubit coupler biases $co,ij ^'^'^ l/^l -compensator 
biases ,- were reset to provide the desired Kij and hi, 
respectively. As a consequence of the changing interqubit 
coupler susceptibility, the i-tuner biases also had to be 



updated to hold the qubit inductance constant, as de- 
scribed in Rcf. 28. Finally, the effects of junction asym- 
metry in the CJJ interqubit couplers and \lq \ compen- 
sators led to a small offset of each qubit 's degeneracy 
point (5$°, as shown in Ref. 29. These offsets were in- 
dependently calibrated for each qubit prior to running 
the processor. To compensate, the PMM controlling the 
qubit flux offset was updated. 




FIG. 7: (Color online) Depiction of the measurement wave 
form sequence. Only time-dependent biases have been shown. 
All qubits are subjected to common $ccjj(i) a-nd Ig{t). All 
QFFs are subjected to a common $fatch(^)- Example dc- 
SQUID flux and current iro{t) wave forms shown only 
for device connected to go. 

With the details concerning the PMM aside, the analog 
control signal sequence for implementing the AQO algo- 
rithm introduced in Sect. II was quite simple, as depicted 
in Fig. 7. During the annealing phase there were only two 
time-dependent signals applied to the entire unit cell: the 
CCJJ control signal $^j,jj(t) and the IJ^I -compensation 
signal Ig(t). For simplicity, we chose to use a relatively 
slow 100-/is-long linear ramp of ^^^^^{t) through the re- 
gion relevant for qubit operation between $™jj/<I'o = 
-0.590 and ^^^-^/^q = -0.650, where the superscripts 
m and b denote biases for which the CCJJ rf-SQUID 
potential energy is deep in the monostable and bistable 
regime, respectively. The aforementioned ramp time was 
sufficiently long to avoid the loss of ground state proba- 
bility to Landau-Zener transitions at anticrossings'^^'^'^'*''. 
Experimental data and numerical simulations support- 
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ing this latter statement will be discussed in Sect. V. 
The time-dependent wave form Ig{t) = a\l^{t)\, with 
a = 37.6, was constructed by sampling the modeled 



ccjj(0)|; shown as the solid curve in Fig. 4(a), at 



200 points between $™ j and ^'^^^y 

At the end of the annealing phase, the final task that 
involved time-dependent wave forms was that of reading 
out the state of the unit cell. The states of all eight qubits 
were simultaneously latched by raising the QFP tunnel 
barriers in unison using their global bias ^fatchC^)^^'^*^- 
Thereafter, the dc-SQUID current biases and flux biases 
were modulated to serially read out the state of each 
QFP. In Fig. 7, only those current and flux bias wave 
forms for operating the readout connected to go have 
been depicted. The entire wave form sequence, including 
programming of the 992 PMM elements, a 1 ms cooling 
period following PMM programming, and 1000 repeti- 
tions of the anneal and read sequence depicted in Fig. 7 
was completed in w 2.5 s. 

The key experimental parameter that was varied in 
this study was the temperature of the processor T. We 
calibrated T in situ using macroscopic resonant tunneling 
(MRT) measurements versus $1^ at $|?^jj/$o = -0.6238 
using a dwell time of 1 ms on qubits qo and (74^^ . At this 
particular CCJJ bias, the qubits reached thermal equi- 
librium within the specified dwell time and the resultant 
population statistics versus <J>^ could be fit to the form 



-Pmrt — 2 



1 + tanh 



2\IP\ ($^-$°) 
2kBT 



where \l^ 



1.26 ± 0.02 fj,A was obtained from the inde- 
pendent calibration shown in Fig. 4(a) and was the 
independently calibrated qubit flux offset. This then left 
only one free parameter, namely T. 

The unit cell was tested using the set of 6400 prob- 
lem instances at three values of T. Each instance was 
run 1000 times and the probabilities of observing each 
of the 256 possible \s) were recorded. Given the small 
size of the problems it was trivial to independently de- 
termine |sgs)- Histograms of the probability of observing 
the ground state Pgs are shown in Fig. 8. The output of 
numerical simulations shown in these plots will be dis- 
cussed in Sect. V. Here, the vertical axes correspond to 
the number of problem instances for which the hardware 
returned Pgs within the bounds of each bin of the his- 
togram. The data taken at T = 20 mK, see Fig. 8(a), 
show an obvious spike in the highest bin, 0.95 < Pgs < 1, 
and a tail that appears to be exponentially suppressed at 
lower Pgs. For T = 35 mK, see Fig. 8(b), the spike in the 
highest bin appears diminished and a broad hump in the 
vicinity of Pgs ~ 0.6 is observed. Finally, for T = 50 mK, 
see Fig. 8(c), the high Pgs spike has been flattened and 
the broad hump has shifted to slightly lower Pgs 0.5. 
Nonetheless, for all three values of T, Pgs > 0.1 for all 
instances. 

To gain further insight into the source of the structure 
of the histograms shown in Fig. 8, we selected an exem- 
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FIG. 8: (Color online) Histograms of the experimentally ob- 
served and the simulated probability of finding the system 
in the ground state Pga from a test involving 6400 unique 
random ISG instances at three different temperatures: (a) 
T = 20 mK, (b) T = 35 mK, and (c) T = 50 mK. 



plary problem instance from each of the three features 
in the T = 35 mK results: Example A from the spike 
(Pgs ^ 1), Example B from the broad hump (Pgs ^ 0.6), 
and Example C from the tail (Pgs ^ 0.4). The in- 
stance settings hi and Kij for these three examples are 
given in Appendix B. We then ran each of these ex- 
ample instances 1000 times at a series of temperatures 
20 mK < T < 100 mK. A summary of the observed pop- 
ulation statistics P is shown in Fig. 9. The output of 
numerical simulations shown in these plots will be dis- 
cussed in Sect. V. For Examples A, B, and C, a total 
of 28, 13, and 38 states were observed above the read- 
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FIG. 9: (Color online) Experimental (symbols) and simu- 
lated (curves) probabilities of finding the system in low energy 
states, (a) Example A. (b) Example B. (c) Example C. Solid 
(dashed) curves have been color-coded to match filled (hol- 
low) symbols. 



The results have been summarized in Fig. 10. This plot 
reveals that the ground state |sgs) = |101), |200), and 
1 97) for Example A, B, and C, respectively. Likewise, the 
first excited state(s) can be identified as |s) = |69), |232), 
and the degenerate pair (|54),|101)) for Example A, B, 
and C, respectively. From a comparison of Figs. 9 and 
10 one can make two important observations: First, for 
all three examples the order of the observed states, from 
most to least probable, exactly matches the order when 
ranked by ascending E{s). Thus the ground state is the 
most probable state, the first excited state is the second 
most probable, and so forth. The only exception is the 
order of |240) and (|45),|100)) in Example A, though this 
discrepancy is within the noise. Second, there is a pre- 
ferred order for some of the pairs of degenerate excited 
states such as (|54),|101)) and (|190),|117)) for Example 
C. While the preference for one over the other appears 
to be subtle, it is nonetheless very reproducible. 
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FIG. 10: (Color online) Evaluation of the dimensionless ob- 
jective function E{s) [Eq. (1)] for Examples A, B, and C. 
Instance settings are given in Appendix B. Results plotted 
relative to the lowest value iJ(sgs). Spin configurations \s) are 
identified to the right of each level. 

To summarize this section of the paper, we have tested 
an cight-qubit unit cell inside a superconducting proces- 
sor that embodies an AQO algorithm. The unit cell has 
been tested by running 6400 random ISG problem in- 
stances at three temperatures and three representative 
instances at multiple temperatures. With repeated mea- 
surements, the hardware returns a small number of final 
spin configurations, with the ground state always being 
the most probable result. The remaining states that are 
observed correspond to low lying excited states. 



V. SIMULATED PROCESSOR PERFORMANCE 



out resolution threshold at T = 91.7 mK. respectively. 
We have plotted the data for only those states that were 
observed with P > 0.01 at T = 91.7 mK. 

The results shown in Fig. 9 indicate that a small num- 
ber of |s) are returned for any given problem instance 
and that the probability of observing those states evolves 
smoothly with T. To assist in determining why particu- 
lar states were observed, we have evaluated the objective 
function E{s), given by Eq. (1), for Examples A, B, and 
C using the values of hi and Kij provided in Appendix B. 



The objective of this section of the paper is to present a 
physical picture that naturally explains the final distribu- 
tion of \s) that was obtained from any given ISG problem 
that was posed to the hardware. This picture relies upon 
having independently calibrated all qubit, coupler, and 
I I -compensator parameters that enter into Hamilto- 
nian (6a) and then adding interactions between all eight 
qubits and an environment. It will be demonstrated that 
the inclusion of quantum mechanical relaxation*^ is suf- 
ficient to explain the observed distributions. 
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We have chosen to model this system by neglecting 
all but the two lowest lying states in each CCJJ rf 
SQUID (qubit approximation). Furthermore, we have 
also assumed that the qubit parameters | J|'($^(.jj) | and 
Ag($^^jj) are the same for all qubits and are given by 
the solid curves in Fig. 4. This latter assumption would 
not be justified without having first leveraged all of the 
features of the CCJJ rf SQUID so as to synchronize their 
qubit parameters to a reasonable level. The low energy 
physics of the ideal closed system will then be dictated 
by Hamiltonian (6a). 

Instantaneous eigenspectra as a function of annealing 
parameter ^ccjj ^^^^ example problem instances A, B, and 
C, as generated using Hamiltonian (6a), are shown in 
Fig. 11. Note the negative sign in front of the indepen- 
dent variable in these plots, which has been inserted such 
that annealing progresses from left to right. The eigcnen- 
ergies E have been plotted relative to the ground state 
energy Eq such that the ground state is coincident with 
the horizontal axis. The lowest 16 levels have been plot- 
ted for Examples A and B and the lowest 32 levels for 
Example C. In all cases the eigenspectra have been trun- 
cated above the second excited state of an isolated CCJJ 
rf SQUID, here denoted as i?rfs, as Hamiltonian (6a) is 
not the correct representation of the physics of the closed 
system beyond this energy. This latter point is of signifi- 
cant practical interest and will receive more attention in 
Sect. VI. The initial states have been identified on the 
left of each plot. Here, the ground state is denoted as 
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and is a product of the isolated qubit ground states. 
There is an octet of first excited states, denoted as 
|ni,m), which are product states of all qubits in their 
ground state except a single qubit m in its excited state, 
(l~l)m ~ l+l)m) /v^- Analogously, there are 64 second 
excited states, denoted as IH2. in which all but two 
qubits are found in their ground state. In all cases, the 
excited states drop in energy with increasing annealing 
parameter until they reach a minimum in the vicinity 
of $^(.jj/$o ^ —0.615. Each of the minima correspond 
to a point in the annealing process at which the ground 
state, initially |Ho), sheds a localized spin configuration 
|s) from the superposition, thus eventually localizing in 

the solution state |sgs). Likewise, the excited states be- population statistics appear to freeze during annealing, 
come localized in configurations |s) 7^ |sgs) as they are 
ejected from the ground state, and then increase in en- 
ergy to the right of the minima. In a large-scale QSG, 
this process would give rise to a second-order phase tran- 
sition between a ground state possessing a high degree of 
symmetry (a superposition of all |s)) to the left of the 
minima and a ground state with broken symmetry (lo- 
calized in |sgs)) to the right of the minima. Those |s) 
that were observed with P>0.01atr = 91.7 mK, as in- 
dicated in Fig. 9, have been labeled on the right of each 
corresponding plot in Fig. 11. 



FIG. 11: (Color online) Calculated low energy eigenspectra 
versus — (a) Example A. (b) Example B. (c) Example 
C. Energies have been plotted relative to the ground state 
energy Eq. Energies have been truncated above that of the 
second excited state of the CCJJ rf SQUID, _Erfs — i?o- Initial 
states are labeled as lllo), |ni_m) and |n2,m,n) on the left. 
Final states are labeled as \s) on the right. Vertical dashed 
lines roughly delineate the value of — $ccjj beyond which the 



Each plot in Fig. 11 has been marked with a vertical 
dashed line at $|?(,jj/$o = —0.622 to roughly delineate 
where in the anneahng process the population statistics 
appear to freeze. As a crude approximation, one can take 
the energies of the excited states along this vertical line 
and calculate Boltzmann occupation factors that roughly 
agree with the results presented in Fig. 9. However, such 
a simplistic picture cannot explain the very reproducible 
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disparity between the observed probabihties of excited 
states that are degenerate along the vertical line. 

To begin the discussion of our dynamical model, it is 
prudent to briefly address what is not the source of ei- 
ther the structure of the Pgg histograms in Fig. 8 or the 
nonzero probability of observing excited states in Fig. 9. 
We argue that none of these observations can be ex- 
plained by a violation of the adiabatic condition. First, 
we have studied all 800 progenitor random ISG problems 
by numerically solving a dynamical model the system 
described by Hamiltonian (6a) in the absence of environ- 
mental noise and subject to the wave forms depicted in 
Fig. 7. These simulations revealed that the 100 /xs CCJJ 
bias ramp time was sufRciently long as to deem the evo- 
lution of the closed system adiabatic, as they always re- 
turned Pgs = 1 to within numerical error. Second, we 
have experimentally studied Examples A, B, and C at 
multiple CCJJ bias ramp times from 10 /is to 1 ms and 
observed relatively little change in the outcome. Band- 
width limits on our apparatus prevented us from studying 
the system subject to shorter ramp times. 

We now turn to a description of the open system model 
that we have investigated. We wish to first deal with 
the physics of relaxation and then discuss the effects of 
dephasing thereafter. Furthermore, we must carefully 
distinguish between AQO, in which the evolution of the 
closed system can be described as being adiabatic, and 
QA, in which the coupled qubit system exchanges energy 
with its environment. 

It has been argued that flux noise is the principal 
soiurce of decoherence for superconducting qubits of the 
type included in our circuits^^. For solid state qubits, 
where materials defects proximate to the qubit body are 
most likely the intermediaries by which a thermal envi- 
ronment disturbs the state of a qubit, this is best cap- 
tured by coupling each qubit to its own independent local 
environment. The Hamiltonian for this open system can 
be written as 

H{t) = Ho{t) + 2 E 2 Wl Q^iM''' , (10) 

i=0 

where 'Ho{t) is given by Hamiltonian (6a) and Qi(t) is an 
operator that acts upon the environment seen by qubit i 
that has units of flux. Note the factor of \IP {t) \ preced- 
ing each Qi{t), which captures the fact that as annealing 
progresses that the coupling between qubit i and environ- 
mental flux noise will change. The noise spectral density 
for the environment surrounding qubit i as a function of 
angular frequency w can be expressed as^^ 

S.,{oj) = (2 \IP {t)\f j dt'e'^'' (Q.(t')Q.(O)) . (11) 

For Markovian noise*^, the fluctuation-dissipation theo- 
rem allows one to write 



Si{uj) = 



2lui{L^[u)} {2\iv{t) 




1 — exp 1 




' ksT ) 





(12) 



where Lq{uj) can be thought of as an inductance per unit 
bandwidth that is imparted to qubit i by its local en- 
vironment. At a more fundamental level, Lq must be 
proportional to the environment's complex magnetic sus- 
ceptibility summed over all modes. 

While there is general consensus that the flux noise 
power spectral density in superconducting flux qubits is 
1// in nature at very low frequency*®, relatively little 
is known about the form of S'i(w) that is appropriate at 
the high energies required to induce transitions in these 
qubits. In lieu of a more detailed picture of Si(u;), we 
have chosen to follow a simple but physically reasonable 
approach by taking the Ohmic approximation in which 
lm{Lo(uj)} is treated as a constant. This simple quantum 
mechanical model of an environment is sufficient for our 
purposes as it provides an asymmetric Si{io) with respect 
to Lo. Consequently, the fiuctuation-dissipation theorem 
ensures that there are thermally occupied modes of the 
environment that can fuel excitations of the eight-qubit 
system, thus compromising Pgg. 

To model the dynamics of the open system embod- 
ied by Hamiltonian (10), we have constructed a nu- 
merical model based upon a generalized Bloch-Redficld 
formalism*^'^°. In this model, one integrates a set of cou- 
pled first order differential equations for the elements of 
the eight-qubit system's density matrix p{t) , as given by 

+ Mmnki) Pki , (13a) 

k,l 

Mmnki = -5nk {l\ df \m) - 5ml {n\ dt \k) , (13b) 

where Rmnki is the so-called Redficld tensor^° which is 
constructed from matrix elements of the fiux operators 
between the eight-qubit states |m) and the Si{ijj) of the 
individual qubits. These matrix elements, which vary 
proportional to the number of qubit flips between a pair 
of states, are what facilitates thermalization of the eight- 
qubit system. In order to make the numerical work man- 
ageable on a personal computer, we kept only the ten 
lowest lying states. We have verified that the output of 
this model in the limit of small Im{Lo} matches that of 
numerically solving the Schrodinger equation. 

Our numerical model of the open system contains only 
one free parameter, namely, Im{Lo}- One can generate 
a crude estimate of this parameter based upon the width 
of a single qubit lowest order MRT peak deep in the 
incoherent regime. From Eq. (9) of Ref. 51, 

= j dujS,{u) . (14) 

For the qubits in the circuit studied in this paper, ex- 
perimental results yielded Gaussian-shaped MRT peaks 
of width W/2\IP\ - 80/i«>o at r = 20 mK. Low fre- 
quency drift measurements, of the type described in 
Ref. 52, revealed low frequency fiux noise with a one- 
sided power spectral density that could be fit to the 
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form I f with A ~ l/x$o/V^Hz- Integrating this 1// 
power spectrum over a nominal ten decades in frequency 
yields a root mean squared flux noise on the order of 
Wlf/2 I I ~ 5 ^$0- Therefore, there was an additional 
{W - Wlf) /2 [/fl ~ 75 ^$0 of integrated flux noise that 
was not captured by the 1// portion of the noise spec- 
trum. Assuming that the remainder of can be 
attributed to a small white flux noise level that is inte- 
grated out to a cutoff /c = ksT/h ^ 460 MHz. we esti- 
mate that the amplitude of the one-sided flux noise spec- 
tral density will be on the order of (5$„ ^ 3.5 w^q/ ^/ H z. 
This white noise level is well below the detection limits 
realized in all low frequency flux noise measurements on 
SQUIDs reported in the literature to date. Experimen- 
tal evidence of an even lower white flux noise level in 
a 3-junction rf SQUID flux qubit has been reported in 
Ref. 53. In terms of an environmental inductance, our 
estimate of (5$„ translates into Im{Lo} = 47yH/Hz. 

The results of numerical simulations of Eq. (13a) with 
the estimated Im{Lo} and the problem instance settings 
for Examples A, B, and C as a function of T are sum- 
marized in Fig. 9. Here, the coloring and format of the 
curves have been chosen to correspond to the color and 
fill of the data symbols in those plots, with hollow (solid) 
points corresponding to dashed (solid) curves. While the 
match between theory and experiment is not perfect, it 
is nonetheless clear that the generalized Bloch-Redfield 
model does provide a reasonable prediction of the popula- 
tions for all three Examples. We believe that most of the 
discrepancies can be explained by slight errors in prob- 
lem specification, as will be discussed below, and there- 
fore did not pursue any fitting of these data by chang- 
ing Im{Lo}. Very importantly, we draw attention to the 
splitting of the population statistics between the pairs of 
degenerate excited states (|54),|101)) and (|190),|117)) 
for Example C. These features arise quite naturally from 
the model and are a result of differences in the transition 
rates between states during the course of annealing, as 
governed by the quantum mechanical matrix elements of 
the flux operators between pairs of initial and final states. 

In addition to the three examples cited above, we have 
simulated all 6400 problem instances using our model and 
constructed Pgs histograms to compare with the data in 
Fig. 8. The results of these simulations are summarized 
in those plots. These simulations did indeed yield the 
expected general features, specifically the loss of counts 
from the highest 0.95 < Pgs < 1 bin to the broad hump 
centered about Pgs ^ 0.6 with increasing T, and an expo- 
nentially suppressed tail at low Pgs. However, the sim- 
ulations also yielded more small-scale structure in the 
histograms than was observed in experiment. 

It was hypothesized that slight errors in embedding 
problem instances in the hardware could have been re- 
sponsible for the absence of fine structure in the exper- 
imental histograms in Fig. 8. The PMM elements re- 
sponsible for nulling small qubit flux offsets are on-chip 
digital-to-analog flux converters whose least significant 
bit (LSB) weights were measured to be iJ'I'^^^ « 38/x$o- 



Therefore, we cannot null qubit degeneracy point off- 
sets <I>° to better than ±5¥f^ /2 using PMM. Further- 
more, the PMM elements controlling the coupler and 
I /I* I -compensator fiux biases, <&co,y and ^f^^i, had LSB 
weights 5$co^^ ~ 5<^Y^ w 60/i$o- This latter quantity 
limits the precision to which Kij oc My and hi cx Mi 
can be specified, respectively, thereby distorting problem 
instances. We have attempted to simulate the effects 
of the finite precision of our PMM elements on all 6400 
problem instances. It was observed that such embedding 
errors only marginally smoothed out the fine structure in 
the simulated histograms. 

A second potential source of error in embedding prob- 
lem instances onto the hardware was that due to imper- 
fect setting of the CCJJ control biases that are used to 
null the effects of Josephson junction variations^®. Such 
imperfections give rise to an apparent flux offset in a 
qubit that varies monotonically with the annealing con- 
trol, meaning that — ^ <i>° (•I'ccjj) in Eq- (6b), as given 
by Eq. (B4c) of Ref. 28. Our single-qubit calibration pro- 
cedures null |$o - $0 I to less than 20 ^$o, 
where the biases and ^f!^,: straddle the region that 
is critical for annealing, as depicted in Fig. 7. Thus, while 
such residual imperfections will inevitably be present, we 
nonetheless consider their impact to be subdominant to 
the errors introduced by the flnite precision of the PMM. 

Our multiqubit generalized Bloch-Redfleld model in- 
vokes Ohmic noise spectral densities that extend to high 
energy, thus providing a mechanism by which each indi- 
vidual qubit can exchange energy with its environment. 
Within the context of gate model quantum computation 
(GMQC), high frequency noise of this type gives rise to 
the timcscalc Ti over which a qubit prepared in its ex- 
cited state would relax towards thermal equilibrium. In 
many respects, QA is compromised by the reverse process 
in which a qubit prepared in its ground state is excited by 
the environment over a timescale e^i/'^^'^Ti, as dictated 
by the Einstein relation between upward and downward 

transition rates, where Eq = \J ^"q ^"q excita- 
tion energy of a qubit for arbitrary bias conditions. This 
implies that QA is somewhat more robust against the 
effects of thermalization than GMQC, although our nu- 
merical simulations clearly indicate that the ground state 
probability Pgs is reduced by these Ti-related processes. 
Nonetheless, theoretical studies such as Ref. 50 have con- 
cluded that the computation time need not be severely 
compromised by the presence of such an environment. 

One of the key oversights of our generalized Bloch- 
Rcdfield model is the role of non-Ohmic low frequency 
flux noise, particularly noise spectral densities of the form 
S'i(w) cx l/o;'' (77 > 0) that are commonly referred to as 
1// noise. Within the context of GMQC, low frequency 
flux noise adiabatically manipulates qubit biases, thus 
leading to dephasing over a timescale frequently denoted 
as Tip . From the perspective of Q A, it is more convenient 
to envision low frequency flux noise as being responsi- 
ble for the generation of random flux offsets within each 
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qubit that change with every repetition of an experiment. 
These flux offsets will possess Gaussian statistics whose 
distribution has a width given by the integrated low fre- 
quency flux noise Wlf/2 j^'^- According to Eq. (6b), 
random flux offsets will distort the individual hi by al- 
tering , thus shifting the intended Hamiltonian off tar- 
get. If the random flux offsets are small, one can still 
expect the hardware to return |sgs) with high Pgg, but 
that probability will differ from predictions made from 
simulating the system in the absence of low frequency 
flux noise. If the random flux offsets are comparable to 
the applied local flux bias on a given qubit at the criti- 
cal point during annealing when the population statistics 
freeze, Wlf/2 \lq\^ ^ hiM afm then the hardware 

could end up 'solving' the wrong problem. Thus, low fre- 
quency noise ultimately impacts the precision to which 
one can specify a problem instance. This issue is cur- 
rently being investigated at a theoretical level. 

In terms of a heuristic picture, it appears that QA, 
when run slowly with respect to the adiabatic limit, pro- 
gresses as follows: First the system is initialized at an 
annealing parameter ^ccjj(^ = 0) = <I>™jj where the gap 
between ground and first excited state is much greater 
than ksT. The large gap ensures that the ground state 
is occupied with near certainty. Next, as ^ccjj is slowly 
lowered the system thermalizes to the environment. The 
environment strives to bring the qubit system to a Boltz- 
mann distribution, though the ultimate achievement of 
that outcome can be hindered by differing transition rates 
between states. The system continues to thermalize un- 
til the qubit dynamics are significantly impeded by the 
reduction of Ag. Thereafter, despite the fact that the 
eigencnergies continue to diverge due to the growth of 
qubit persistent currents, the state of the processor will 
no longer evolve. Thus the population statistics become 
frozen at the levels achieved at some higher value of 

'i>?cjj(t)>1>?cjj(i/) = ejr 

To summarize this section of the paper, we have simu- 
lated the performance of the coupled eight-qubit system 
using a generalized Bloch-Redfield model to capture the 
physics of relaxation processes that couple the individual 
qubits to an environment in thermal equilibrium. The 
amplitude of the white noise spectral density sensed by 
each qubit was estimated from the widths of macroscopic 
resonant tunneling peaks. Comparison of the output of 
the model to the temperature-dependent probability dis- 
tributions for three example problems indicated that the 
model captures the broad features seen in the data. The 
success histograms for all 6400 problem instances run on 
the hardware were roughly reproduced. 



VI. DISCUSSION AND OPEN QUESTIONS 

The experimental results in Sect. IV and the modeling 
in Sect. V present a compelling physical picture of how 
the state of our prototype processor evolves when run 
slowly in comparison to the adiabatic limit. These results 



have prompted some interesting conclusions and a series 
of open questions that we will briefly address below. 

Our study has demonstrated that the most proba- 
ble result from the hardware, when run slowly, is in- 
variably the ground state. The remaining probability 
is distributed between low lying states that are roughly 
within an energy window O^ksT) above the energy of 
the ground state at the point in the annealing process 
where the population statistics freeze. Depending upon 
the nature of the optimization problem, such low lying 
states may still constitute acceptable solutions. This ob- 
servation suggests two possible modes of running a QA 
processor in practice: 

• Run a given optimization problem a statistically 
large number of times. Take all \s) that are 
observed with probability greater than the read- 
out resolution threshold and calculate E{s) using 
Eq. (1). Take whatever \s) provides the lowest E 
as the solution. 

• Run a given optimization problem once. Take the 
output |s) and calculate E{s) using Eq. (1). If E is 
less than some user-defined threshold, then accept 
the solution. If not, then iterate. 

Note that had our processor been operated in the first 
mode, then it would have returned the correct answer to 
all 6400 problem instances. 

We now turn to the open questions that have been 
motivated by this work: 

What happens when the CCJJ bias ramp time is re- 
duced? Given the limited bandwidth of our external bias 
lines, we were unable to elicit a substantial change in pop- 
ulation statistics. Thus, we were restricted to studying 
the regime in which the thermalization times were much 
shorter than the annealing time. Consequently, we can 
only claim to have demonstrated QA, not AQO. Efforts 
are underway to build a new apparatus that will allow us 
to probe the regime in which thermalization times exceed 
the annealing time in an eight-qubit processor. 

How does the performance scale with problem size? 
The eight-qubit unit cell is too small to be used to address 
issues concerning scaling. The intention of this study was 
to provide a basic demonstration of what we believe to be 
a complete set of essential ingredients needed for build- 
ing a scalable QA processor acting in concert. We will 
reserve a discussion of the operation of larger portions of 
a complete 128-qubit chip to a future publication. 

How does one implement error correction in QA ? The 
work presented in this article is in much the same spirit as 
that of Ref. 55 - operate a small scale device that 'looks' 
like a basic quantum information processor and run it 
using the simplest algorithms known. We have made no 
attempts beyond statistical sampling to implement any 
form of error correction in these particular experiments. 
This is an active area of research at the moment. 

To what precision must the qubit parameters 1 1^ | and 
Aq be synchronized as a function of^^^jj{t)? Implicit in 
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many models of AQO is the assumption that all of the 
qubits are identical throughout the annealing process. 
An examination of the data in Fig. 4 reveals that while 
is very uniform as a function of annealing parame- 
ter $ccjjj there are clear discrepancies between values of 
Ag at the same bias While we have demonstrated 

that our modest eight-qubit unit cell is capable of solving 
small-scale ISG problem instances, it is neither clear to 
what degree the asynchronization of A, (t) may have im- 
peded its performance nor what the implications of such 
asynchronization are for large-scale systems. 

How do the higher excited states of the CCJJ rf 
SQUIDs affect performance? To our knowledge, this is- 
sue has not been addressed in the literature. We strongly 
advise against literal interpretations of many-qubit eigcn- 
spectra to energies that exceed that of the second excited 
state of an rf-SQUID qubit. This has consequences for 
studies regarding the ultimate utility of QA. Higher ex- 
cited states provide fast interwell relaxation mechanisms 
that could limit failures during a computation. 

VII. CONCLUSIONS 

An eight-qubit unit cell that is part of a supercon- 
ducting optimization processor has been experimentally 
investigated. This processor makes use of several scalable 
elements, including on-chip programmable flux sources, 
XY-addressablc high fidelity readouts, and the use of a 
limited number of global analog control lines to provide 
a variety of time-dependent control signals to multiple 
devices. The processor was tested using a large set of 
randomly generated Ising spin glass problem instances. 
The experimental results were shown to be consistent 
with the predictions of a quantum mechanical theory in 
which the individual qubits are coupled to a thermal en- 
vironment. 

We thank J. Hihon, F. Brito, K. Pudenz, S. Han, A. 
Kleinsasser, and G. Kcrbcr for useful discussions. 



Appendix A: Glossary of Abbreviations 

AFM antifcrromagnctic 

AQO adiabatic quantum optimization 

CCJJ compound-compound Josephson jimction 

CJJ compound Josephson junction 

GMQC gate model quantum computation 

ICO internal coupler 

IPC persistent current compensator 

ISG Ising spin glass 

LSB least significant bit 

LT inductance tuner 

MKT macroscopic resonant tunneling 

PMM programmable magnetic memory 

QA quantum annealing 

QFP quantum flux parametron 

QSG quantum Ising spin glass 

RO readout 

XCO external coupler 

IQLZ single-qubit Landau-Zener 

2QLZ two-qubit Landau-Zener 



Appendix B: Example Problem Instance Settings 



Parameter 


Example A 


Example B 


Example C 


ho 


-1 


-5/7 


-3/7 


hi 


4/7 


6/7 


2/7 


h2 


5/7 


-5/7 


6/7 




4/7 


-1 


3/7 


hi 


-4/7 


6/7 


-5/7 


/l5 


1/7 


-3/7 


2/7 


he 


-5/7 


-3/7 





h7 


4/7 


-3/7 


-2/7 


Ko4 


2/7 


-6/7 


-5/7 




4/7 


-6/7 


1/7 


K24 


-1/7 










-1/7 


6/7 


-3/7 


Ko5 


4/7 


6/7 


5/7 


Ki5 


-2/7 


1/7 


5/7 


K25 


-1/7 


3/7 


-1/7 


Ks5 


6/7 


-1/7 


-4/7 


K06 


1/7 


1 


-6/7 


Kie 


1/7 


3/7 


5/7 


K26 


4/7 


-5/7 


-5/7 


if 36 


-5/7 


2/7 


-1/7 


Km 


6/7 


1/7 


1/7 


Kn 


2/7 


4/7 


-3/7 


K27 


4/7 


6/7 


-6/7 


K37 


3/7 
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